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ABSTRACT 

We report the discovery of three proplyd-hke structures in the giant HII region NGC 3603. The 
emission nebulae are clearly resolved in narrow-band and broad-band HST/WFPC2 observations in the 
optical and broad-band VLT/ISAAC observations in the near-infrared. All three nebulae are tadpole 
shaped, with the bright ionization front at the head facing the central cluster and a fainter ionization 
front around the tail pointing away from the cluster. Typical sizes are 6,000 A.U. x 20,000 A.U. The 
nebulae share the overall morphology of the proplyds ( "PROto PLanetarY DiskS" ) in Orion, but are 
20 to 30 times larger in size. Additional faint filaments located between the nebulae and the central 
ionizing cluster can be interpreted as bow shocks resulting from the interaction of the fast winds from 
the high-mass stars in the cluster with the evaporation flow from the proplyds. 

Low-resolution spectra of the brightest nebula, which is at a projected separation of 1.3 pc from the 
cluster, reveal that it has the spectral excitation characteristics of an Ultra Compact HII region with 
electron densities well in excess of 10^ cm~'^. The near-infrared data reveal a point-source superimposed 
on the ionization front. 

The striking similarity of the tadpole shaped emission nebulae in NGC 3603 to the proplyds in Orion 
suggests that the physical structure of both types of objects might be the same. We present 2D radiation 
hydrodynamical simulations of an externally illuminated star-disk-envelope system, which was still in its 
main accretion phase when first exposed to ionizing radiation from the central cluster. The simulations 
reproduce the overall morphology of the proplyds in NGC 3603 very well, but also indicate that mass- loss 
rates of up to 10~^ Mq yr~^ are required in order to explain the size of the proplyds. 

Due to these high mass-loss rates, the proplyds in NGC 3603 should only survive wlO^ yr. Despite this 
short survival time, we detect three proplyds. This indicates that circumstellar disks must be common 
around young stars in NGC 3603 and that these particular proplyds have only recently been exposed to 
their present harsh UV environment. 

Subject headings: circumstellar matter - stars: formation - stars: pre-main sequence - open clusters 
and associations: individual (NGC 3603) - ISM: individual (NGC 3603). 



1. INTRODUCTION 

HST/WFPC2 observations of the Orion Nebula (M42) 
revealed a large variety of dark silhouette disks (O'Dell 
& Wong 1996; McCaughrean & O'DeU 1996) and par- 
tially ionized circumstellar clouds (O'Dell et al. 1993). 
Many of the circumstellar clouds, which had first been de- 
tected from the ground by Laques & Vidal (1979), have a 
cometary shape with the tails pointing away from the 07V 



star ©1 Ori C and the 09.5V star 6^ Ori A, the bright- 
est and most massive members of the Trapezium cluster. 
The partially ionized circumstellar clouds with cometary 
shape were identified as protoplanetary disks ( "proplyds" ) 
around young stars, which are ionized from the outside 
(Churchwell et al. 1987; O'DeU et al. 1993). 

Many proplyds are ionization bounded, which indicates 
that all EUV photons (hi' > 13 eV) get absorbed in the 
ionization front engulfing the protostar and its circum- 
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stellar disk (O'Dell 1998). FUV photons (13 eV > hv > 
6eV), however, are able to penetrate the ionization front. 
They heat up the inside of the proplyd envelope and lead 
to the dissociation of molecules in the outer layers of the 
circumstellar disk (Johnstone et al. 1998). The resulting 
evaporation flow provides a steady supply of neutral atoms 
to the ionization front and leads to the development of a 
cometary tail (McCuUough et al. 1995; Storzer & HoUen- 
bach 1999). 

Because of their larger size and the ionized envelope, 
proplyds can be spotted more easily than circumstellar 
disks themselves. Consequently, Stecklum et al. (1998) 
proposed to utilize proplyds as tracers for circumstellar 
disks in distant star forming regions. Systematic search 
efforts for proplyds in HII regions around young clusters 
with WFPC2 did not yield any new detections (Stapelfeldt 
et al. 1997; Bally et al. 1998a). Until recently, only one 
other proplyd had been found. It is located in the vicinity 
of the 07V star Herschel 36 in the Lagoon Nebula (M16, 
Stecklum et al. 1998). 

NGC 3603 is located in the Carina spiral arm at a dis- 
tance of 6kpc (De Pree et al. 1999 and references therein). 
With a bolometric luminosity Lboi > 10^ Lq, NGC 3603 is 
100 times more luminous than the Orion Nebula and has 
about 10% of the luminosity of 30 Doradus in the Large 
Magellanic Cloud (LMC). It is the only Galactic giant HII 
region whose massive central ionizing cluster can also be 
studied at optical wavelengths. The initial mass function 
of the cluster follows a Salpeter type power law with in- 
dex r=— 1.70 for masses > 25 M© and r=— 0.73 for masses 
less than 25 M0 (Eisenhauer et al. 1998), extending from 
Wolf-Rayet stars and 03V stars with masses up to 120 M© 
(Drissen et al. 1995) down to stars of at least 1 M0 (Eisen- 
hauer et al. 1998). The total cluster mass is > 4,000 M©. 

To the south of the cluster is a giant molecular cloud. 
Ionizing radiation and fast stellar winds from the starburst 
cluster are excavating large gaseous pillars. Located about 
20" to the north of the cluster center is the blue supergiant 
Sher 25. This supergiant is unique because its circumstel- 
lar ring and bipolar outflows form an hourglass structure 
similar to that of SN1987A (Brandner et al. 1997a, 1997b). 

As part of a follow-up study on the hourglass struc- 
ture around Sher 25 we observed NGC 3603 with 
HST/WFPC2. In this paper, we report the serendipitous 
discovery of three proplyd-like structures in NGC 3603 
based on HST/WFPC2 and VLT/ISAAC observation and 

perform a first analysis of their physical properties. 
2. OBSERVATIONS AND DATA REDUCTION 

2.1. HST/WFPC2 observations 

On March 5, 1999 we obtained deep narrow-band Ha 
(F656N, 2 X 500s) and [Nil] (F658N, 2 x 600s) observations 
of the giant HII region NGC 3603. The Planetary Camera 
(PC) chip was centered on the bipolar outflow structure 
around the blue supergiant Sher 25. The three Wide Field 
Camera (WF) chips covered the central cluster and the 
HII region to the south of Sher 25. 

In addition, we retrieved and analyzed archival HST 
data, which had originally been obtained in July 1997 (PI 
Drissen). The PC was centered on the cluster, and the 
three WF chips covered the area north-west of the clus- 
ter. Using IRAFF], we combined individual short exposure 

^IRAF is distributed by the National Optical Astronomy Observatories (NOAO) 



in F547M (8x30s), F675W (8x20s), and F814W (8x20s) 
to produce images with effective exposure times of 240s, 
160s, and 160s, respectively. 

The surface brightness of the proplyds was measured us- 
ing aperture photometry with an aperture radius of 0.5". 
The photometric calibration was carried out following the 
steps outlined in the HST Data Handbook Version 3. No 
attempt was made to correct for the contribution of the 
[Nil] lines to the Ha F656N filter, or the contribution of 
the Ha line to the [Nil] F658N filter. The spectrum of 
Proplyd 1 (see below) indicates that the underlying con- 
tinuum emission from the proplyd is negligible. Applying 
equation (3) from O'Dell (1998) yields that the contribu- 
tion of the [Nil] lines to the total flux observed in the Ha 
F656N filter is at most 3.5%. 

2.2. Preparatory ground-based observations 

A first set of deep ground-based broad and narrow-band 
images of NGC 3603 was obtained on April 22, 1991 with 
the ESO New Technology Telescope and the ESO Multi- 
Mode Instrument (EMMI) . These data were used to iden- 
tify a number of compact emission nebula in the vicinity 
of the central cluster. 

On April 2, 1994 we tried to resolve the inner struc- 
ture of the compact nebulae using the ESO Adaptive Op- 
tics system ADONIS. This attempt failed due to the lack 
of sufficiently bright stars suitable for wavefront sensing 
within 20" of any of the proplyds. 

A low-dispersion spectrum of Proplyd 1 was obtained 
on February 3, 1997 with the ESO/MPI 2.2m telescope 
and the ESO Faint Object Spectrograph 2 (EF0SC2) at 
La Silla, Chile. The slit width was 1.5". The spectrum has 
a sampling of 0.2 nm pixel" ^, a spectral resolution around 
450 kms~^, and covers the wavelength range from 517 nm 
to 928 nm. It was wavelength and flux calibrated using 
IRAF. We did not try to correct for fringes, which become 
noticeable redward of 750 nm. 

2.3. VLT/ISAAC observations 

As part of a study of the low-mass stellar content of the 
starburst cluster (see Brandl et al. 1999), NGC 3603 was 
observed with the ESO Very Large Telescope (VLT) Unit 
Telescope 1 (UTl, now officially named "ANTU") during 
the nights of April 4-6 and 9, 1999. The observations 
were carried out in service mode and used the Infrared 
Spectrograph And Array Camera (ISAAC, see Moorwood 
et al. 1998). Deep near-infrared observations of NGC 3603 
were obtained with effective exposure times in Jj, H, and 
Ks of 2,230s, 2,710s, and 2,890s, respectively. The see- 
ing (FWHM) on the co-added frames was of the order of 
0.35" to 0.40". Dithering between individual exposures in- 
creased the field of view from its nominal value of 2.5' x 2.5' 
to 3.5' X 3.5'. The data were flux calibrated based on ob- 
servations of faint near-infrared standard stars from the 
hsts by Hunt et al. (1998) and Persson et al. (1998). 

More details on the data reduction and analysis can be 
found in Brandl et al. (1999). 



3. PHYSICAL PROPERTIES OF PROPLYDS 
3.1. Morphology & Size 
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Fig. 1. — WFPC2 observations of NGC 3603. North is up and east is to the left. The upper part of the image consists 
of the archive data with the foUowing color coding: F547M (blue), F675W (green), F814W (red). Overlaid are our new 
WFPC2 data with the F656N data in the red channel, the average of F656N and F658N in the green channel, and F658N 
in the blue channel. The location of the three proplyd-like emission nebulae is indicated. The insert at the lower right is 
a combination of WFPC2 F656N (blue) and F658N (green) and VLT/ISAAC K^ (red) observations. 



The HST/WFPC2 observations are presented in Fig- 
ure IT]. The figure shows an overlay of two composite 
color images. The upper part of the image consists of 
the archive data with the following color coding: F547M 
(blue), F675W (green), F814W (red). Overlaid are our 
new WFPC2 data with the F656N data in the red chan- 
nel, the average of F656N and F658N in the green channel, 
and F658N in the blue channel. The locations of the pro- 
plyds are marked by small boxes, and enlargements of the 
boxes are shown in the upper part of Figure n. Proplyd 
3 has only been observed in intermediate and broad-band 
filters and thus stands out less clearly against the under- 
lying background when compared to Proplyds 1 and 2. 
The insert at the lower right shows a color composite of 



HST/WFPC2 F656N (blue) and F658N (green) data and 
VLT/ISAAC Kg data (red). 

All three proplyds are tadpole shaped and rim bright- 
ened, with the extended tails facing away from the star- 
burst cluster. The portion of the ionized rims pointing 
towards the cluster are brighter than the rims on the op- 
posite side. The central parts of the proplyds are fainter 
than the rims, with a noticeable drop in surface brightness 
between the head and the tail. 

Proplyds 2 and 3 exhibit a largely axisymmetric mor- 
phology, whereas Proplyd 1, which is also the one closest 
to the cluster, has a more complex structure. Unlike the 
convex shape of the heads of the other proplyds, Proplyd 
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1 has a heart-shaped head with a coUimated, outflow-hke 
structure in between. One possible explanation for the 
more complex morphology of Proplyd 1 might be that it 
is actually a superposition of two (or maybe even three) 
individual proplyds or that the photoevaporative flows of 
several disks in a multiple system interact to produce this 
complex single structure. 

At distances of 7.4" and 2.9" from Proplyd 1 and 2, re- 
spectively, faint arc-like Ha emission features are seen on 
the WFPC2 frames. The arcs are located in the direction 
of the cluster, and may be the signatures of bow shocks 
created by the interaction of proplyd winds with the winds 
from the massive stars in the central cluster. 

The proplyd heads have diameters between 1.2" and 
1.7" (7,200 and 10,800 A.U.). The head-to-tail extent of 
the proplyds is between 2.5" and 3.5" (15,000 to 21,000 
A.U.). In Orion the typical diameters of the proplyd heads 
vary from 45 to 355 A.U. (O'DcU 1998), and the proplyd 
head in M8 has a diameter of 1,080 A.U. (Stecklum et 
al. 1998). Thus, the proplyds in NGC 3603 are 20 to 30 
times larger than the largest proplyds in Orion, and 7 to 
10 times larger than the proplyd in M8. It should be noted 
that proplyds with sizes similar to those of the Orion pro- 
plyds would be too small to be resolvable at the distance 
of NGC 3603, where one pixel (0.1") on the wide field 
CCDs of WFPC2 corresponds to 600 A.U. The PC data 
with a finer pixel scale of 0.0456" per pixel (270 A.U.) re- 
veal indeed several faint point sources which appear to be 
brighter in Ha than in [Nil] . A detailed analysis of these 
sources will be subject of a later paper (Grebel et al., in 
prep.). 

In Orion, the size of the proplyds loosely scales with dis- 
tance from the ionizing source in the sense that proplyds 
further away from O^ Ori C are larger (McCuUough et al. 
1995; Johnstone et al. 1998, O'Dell 1998). In NGC 3603, 
there is no such correlation between the size of a proplyd 
and its projected distance from the cluster. If, however, 
the complex structure of Proplyd 1 results from multiple 
photoevaporating disks as discussed above, the size esti- 
mates based on isolated proplyds cannot be applied. Only 
if Proplyd 1 can actually be decomposed into individual, 
isolated proplyds with diameters around 0.9", would there 
be a tendency for increasing proplyd size with increasing 
distance from the cluster. 

Coordinates, distance from the cluster center, and ap- 
proximate size of the proplyds are given in Table n. 

3.2. Surface Brightness of Proplyds 

The WFPC2 observations reveal that only the outer- 
most layer of Proplyd 1 is ionized, whereas the interior re- 
mains neutral. Table ^ gives the Ha flux and the surface 
brightness of the proplyds as measured from the WFPC2 
frames. No H/3 observations were available which would 
have allowed us to determine the extinction towards the 
proplyds based on the Balmer decrement. Literature val- 
ues for the foreground extinction towards NGC 3603 range 
from Av = 4'"'^« to 5™''s (Moffat 1983; Melnick et al. 1989). 
The Ha flux has thus been corrected for an assumed fore- 
ground extinction of Ana = A"^'^^. The values for the sur- 
face brightness have not been corrected for extinction. 

Proplyds in Orion get fainter with increasing distance 
from the ionization source. Proplyd 2 is about a factor 
of 2.8 fainter than Proplyd 1. If the projected separa- 



tion from the cluster center is comparable to the phys- 
ical distance, then Proplyd 2 should receive a factor of 
(72.5"/43.6")2 = 2.8 less UV photons than Proplyd 1. 
The remarkably good agreement between the number of 
infalling UV photons and the brightness of the proplyd 
suggests that both proplyds are ionization bounded and 
receive most of the ionizing UV photons directly from the 
cluster. 

Proplyd 3 was only observed with intermediate 
and broadband filters. Its red colors (V— R=3.26™, 
R— 1=0.67™) are caused by a combination of the absence 
of any strong emission lines in the passband of the F547M 
filter, foreground extinction, and possibly the presence of 
an embedded central continuum source. 

Results from the VLT near-infrared broad-band pho- 
tometry can be found in Table ^. The surface brightness 
of all three proplyds increases from Jg to H to Kg. The 
photometry for Proplyd 1 and 3 has not been corrected 
for the contamination by the nearby point sources. 

3.3. Nearby Point Sources 

The deep near-infrared observations with the VLT re- 
veal a large number of faint, red point sources. Two point 
sources are detected close to the head of Proplyd 1. As 
can be seen in the lower right insert of Figure 0, one of 
the point sources coincides with the location of the ioniza- 
tion front and may be physically associated with Proplyd 
1. Its red near- infrared colors (see Table ||) indicate that 
the source is highly embedded and/or has a strong intrin- 
sic IR excess. This is in agreement with what one would 
expect for a young stellar object (YSO) surrounded by a 
circumstellar disk and an infalling envelope. Comparing 
with theoretical pre-main-sequence evolutionary tracks by 
Palla & Stabler (1993), we note that a 1 Myr old YSO with 
J=15.4mag and located in NGC 3603 should have a mass 
around 3M0 (Eisenhauer et al. 1998; Brandl et al. 1999). 
The effects of additional local extinction and accretion lu- 
minosity would of course alter this value. In any case it 
should be kept in mind that the relatively high density of 
field sources makes the physical association between this 
point source and Proplyd 1 somewhat uncertain. 

The point source close to the head of Proplyd 3 is al- 
ready detected on the broad-band HST/WFPC2 observa- 
tions (see Figure ||). The WFPC2 images show that the 
point source is actually located in front of Proplyd 3 and 
thus very likely not physically associated with it. 

The Jg, H, and Kg magnitudes of the point sources as- 
sociated with Proplyd 1 and 3 can also be found in Table 
H No central point source is detected in Proplyd 2 with a 
limiting magnitude of Kg < 18.0 mag. 

3.4. Optical Spectroscopy 

Figure 2 shows the low-resolution spectrum of Proplyd 
1. The most prominent emission lines are identified. If 
one assumes an electron temperature of 10** K, which is 
quite typical for HII regions, the flux ratio between the 
[Sii] lines at 671. 7nm and at 673.1 nm yields an electron 
density well in excess of 10"* cm^'^. Unfortunately, such 
densities are close to the collisional de-excitation limit of 
the [Sii] doublet, which prevents us form getting a more 
precise estimate. Other density indicators such as the ex- 
tinction corrected Ha surface brightness or the ratio of the 
[Cm] UV doublet would have to be employed in order to 



i 



Brandner et al. 



-14, 



3.0'10 



2.5'10"^^- 



2.0'10'''^- 



E 



-14 



1.5'10 



LO-IO"""^- ^ 



5.0«10"''^- 



+ 



+ 
O 



600 



X 






Paschen Series 



700 



^ J 



5i[nm] 



J I I L 

800 




/w 



900 



W" 



Fig. 2. — Spectrum of Proplyd 1. The most prominent emissison lines have been marked. The rising continuum towards 
longer wavelengths indicates the presence of an embedded continuum source. 



derive a more accurate estimate of the density (see Henney 
& O'Dell 1999). Densities of the order of 10^ cm'^ to lO*' 
cm~^ have been found in the Orion proplyds (Henney & 
O'Dell 1999), and are also typical for Ultra Compact HII 
regions. 

The infrared source in Proplyd 1 is also detected as 
an underlying, heavily reddened continuum source in the 
spectrum. A more detailed analysis of ground based opti- 
cal images and spectra of the proplyds and other compact 
nebulae in NGC 3603 will be presented in Dottori et al. 
(2000). 

3.5. Limit on the radio brightness 

The proplyds in Orion have been detected with the Very 
Large Array as radio continuum sources at 2 cm and 20 cm 
(Churchwell et al. 1987; FeUi et al. 1993; McCuUough et 
al. 1995). 

NGC 3603 has recently been studied with the Australia 
Telescope Compact Array (ATCA) at 3.4 cm radio contin- 
uum and at several recombination emission lines (De Free 
et al. 1999). Down to the 5 a level of 55 mJy beam~^ for 
a beamsize of 7" none of the proplyds are detected in the 
continuum. 

Following the derivation by McCuUough et al. (1995), 
Section 5, the expected ratio between the radio flux den- 
sity due to Bremsstrahlung, and the Ha flux is 



Fradio _ „ /, r, -O.lrT^O.55 



mJy 



iHa ''"'^* photons cm-2 s-1 ' 

where ^ghz is the radio frequency in GHz and T4 is 



the electron temperature of the ionized gas divided by 
10,000K. 

Thus, for the Ha flux values given in Table || and T4=l, 
one would expect 3.4 cm (8.8 GHz) radio flux densities of 
1.6 mJy and 0.55 mJy for proplyds 1 and 2, respectively. 
This is well below the detection limit of the ATCA obser- 
vations by De Free et al. (1999). Internal absorption of Ha 
photons emitted from the far side of the proplyd amounts 
to on average less than 25% to the total flux (McCuUough 
1993) and has thus been neglected in the above estimate. 

These analytical estimates are in good quantitative 
agreement with the results from the numerical simulations 
described in Section 5. 

4. ENVIRONMENT 
4.1. UV radiation field 

NGC 3603, containing over 20 O stars and WR stars, 
creates a much more extreme UV environment than the 
Trapezium system. The central cluster in NGC 3603 has 
a Lyman continuum flux of 10^^ s~^ (Kennicutt 1984; Dris- 
sen et al. 1995), about 100 times the ionizing power of the 
Trapezium system. 

At the same time, however, the proplyds in NGC 3603 
are at larger distances from the ionizing source. The pro- 
jected separations between the proplyds and the cluster 
center range from 1.3 pc to 2.2 pc. The separation be- 
tween O^ Ori C and the Orion proplyds studied by O'Dell 
(1998) varies between 0.01 pc and 0.15 pc. On average, the 
proplyds in NGC 3603 are exposed to a somewhat less in- 
tense EUV (hi/> 13 eV) radiation field than the proplyds 
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in Orion. 

The most massive stars in NGC 3603 are 03V stars and 
WR stars (Moffat et al. 1994; Drissen et al. 1995), which 
are of earUer spectral type than the late O stars in the 
Trapezium system. As a consequence, the spectral char- 
acteristics of the UV field in NGC 3603 are different from 
those of the UV field in Orion. The early 0-type stars in 
NGC 3603, although more luminous than the late 0-type 
stars in the Trapezium system, produce a disproportionate 
smaller rate of FUV (13 eV > hz^ > 6eV) photons com- 
pared to the rate of EUV photons. Table || gives the ratio 
of the FUV to EUV photon rates for blackbodies with 
temperatures between 30,000 K and 45,000 K. Model cal- 
culations of the atmospheres of hot stars made available 
in electronic form by Adalbert Pauldrachg (see also Paul- 
drach et al. 1998) indicate that for effective temperatures 
below 45,000 K the ratio of FUV to EUV photon rates 
is considerably higher than the ratio derived for a black- 
body. The models indicate ratios of 6.9:1 and 2.4:1 for O 
dwarfs with solar metallicity for effective temperatures of 
35,000K and 40,000K, respectively (see Table |). 

Because all the EUV photons get absorbed in the ion- 
ization front, the FUV photon rate determines the heating 
inside the proplyd envelope and thus ultimately the mass- 
loss rate (Johnstone et al. 1998). 

4.2. Winds and mass-loss rates 

Another important constituent of the environment of 
massive stars are fast stellar winds and wind-wind interac- 
tions. In Orion, four of the five proplyds closest to O^ Ori 
C show arc-like features which may be bow shocks result- 
ing from the interaction of the evaporation flow of the 
proplyds with the fast stellar wind from 8^ Ori C (Mc- 
CuUough et al. 1995). The WFPC2 observations in Ha 
of NGC 3603 reveal similar arc-like features in front of 
Proplyd 1 and 2, that might also be bow shocks. 

For a stationary shock, pressure equilibrium exists on 
both sides of the shock: 



P ^ Pciirci)vli = pp,{rp,)v 



:;2 



where p(r) is the density of the wind at a distance r from 
the cluster (cl) or the proplyd (pr), respectively, and v is 
the terminal velocity of the freely expanding wind. 

If radiative cooling is not important, the stagnation 
point of the shock along the line connecting the cluster 
with the proplyd relates to Adv like 

McWci r^] 



MprWpr 



where M is the mass-loss rate of the cluster and the pro- 
plyd, respectively (e.g., Kallrath 1991). 

Mass-loss rates of individual O stars in NGC 3603 are 
of the order of a few 10~^Mq yr~^, and wind velocities are 
of the order of a few 1,000 kms~^ (Pauldrach et al. 1998). 
If we assume that the 20 O and WR stars in NGC 3603 
have average mass-loss rates of 3xlO~^M0 yr~^, and av- 
erage wind velocities of v=2,000kms^^, the resulting rate 
of momentum input Mv from the combined winds would 
be O.O12M0yr-i kms'^ 

^http://www. USUI. uni- muenchen.de/people/adi/adi. html 



The separation between the arc in front of Proplyd 1 and 
the cluster center is 36.2", and the separation between the 
arc and the proplyd is 7.4". For a bow shock, the ratio 
between the product of mass-loss rate times the velocity 
of the wind from the cluster and Proplyd 1 would thus be 
24:1. 

For a wind velocity of the evaporation flow from Proplyd 
1 of, e.g., 25kms~^, the mass-loss rate would then have to 
be 2x 1O~^M0 yr~^. Whereas this is somewhat on the high 
end side of the parameter space (see section on radiation 
hydrodynamical simulations below), the arc-like feature 
between Proplyd 1 and the cluster could indeed be the re- 
sult of a bow shock. It should be kept in mind, however, 
that the mass-loss rates for the WR stars in NGC 3603 are 
highly uncertain, and might be considerably higher than 
our present estimate. 

For Proplyd 2 the separation between the arc and the 
proplyd is 2.9", and the separation between the arc and 
the cluster center is 69.6". Hence the ratio between the 
product of mass-loss rate times the velocity of the wind 
from the cluster and Proplyd 2 is 580:1, which requires 
much less extreme conditions. For example, a mass-loss 
rate of 2xlO~^MQyr~^ and a wind velocity of lOkms"^ 
would balance the wind force from the cluster at a dis- 
tance of 2.9" from Proplyd 2. Such values are in the range 
of mass-loss rates and flow velocities observed for the pro- 
plyds in Orion (e.g., Hcnney & O'DcU 1999; Bally et al. 
1998b). 

5. RADIATION HYDRODYNAMICAL SIMULATIONS 

The proplyds in NGC 3603 are much larger than the 
proplyds in Orion. In order to investigate if the proplyds 
in NGC 3603 can be explained by a similar physical mech- 
anism as the proplyds in Orion, we carried out radiation 
hydrodynamical simulations. 

5.1. Numerical method and initial conditions 

The simulations are based on a 2D radiation hydrody- 
namics code (Yorke & Welz 1996; Richling & Yorke 1997, 
1998; Richling 1998) and include a thin disk with a finite 
scale height. Diffuse EUV and FUV radiation fields are 
treated in the flux-limited diffusion approximation (Lev- 
ermore & Pamraning 1981) as implemented for multiple 
nested grids by Yorke & Kaisig (1995). The present sim- 
ulation utilizes 6 nested grids with a resolution from 7.3 
A.U. to 233 A.U. The coarsest outermost grid covered a 
cylindrical volume of radius 13,500 A.U. and height 27,000 
A.U. 

The radial and vertical density structure of the disk was 
derived from the collapse of a rotating molecular cloud 
core with a mass of 2 Mq . The collapse was followed up 
through SxlO'^yr, at which time 1.14M0 had already been 
accreted by the central protostar, while 0.86 Mq still re- 
mained in the protostellar disk and the infalling envelope. 
Angular momentum transport was considered via an a- 
prescription (Shakura & Sunyaev 1973), as described by 
Yorke & Bodenheimer (1999). The inclusion of the effects 
of angular momentum transport is important, as disk ma- 
terial in the outer region of the disk actually gains angular 
momentum during the evolution, which results in a very 
extended disk. 

The resulting, initial disk had a diameter of 3,400 A.U. 
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Fig. 3. — Time evolution of a proplyd based on 2D radiation hydrodynamic simulations. The plots show the density 
distribution and evaporation flow velocity field at t=2,416yr, t«50,000yr, and t«138,000yr. The disk can be seen as the 
light structure (high density) in the center of each plot. The ionization front is indicated by the dark solid line. Despite 
an initial mass of 0.8 Mq in the disk and infalling envelope, the evaporation flow with a mass-loss rate of 10~^ Mq yr~^ 
results in an almost complete evaporation of the disk within sa 10^ yr. 



This might seem large compared to the typical disk sizes 
observed for the proplyds in Orion. Based on their nu- 
merical models, Johnstone et al. (1998) derive disk sizes 
between 27 A. U. and 175 A. U. for the proplyds in Orion^. 

For disks which do not show any sign of external il- 
lumination by UV photons, typical disk sizes are in the 
range of 200 to 1,000 A.U. (McCaughrean & O'Dell 1996; 
Padgett et al. 1999). It suggests that the extreme UV ra- 
diation field in the Trapezium cluster is evaporating the 
disks away at a rapid pace. This is in agreement with the 
small remnant disk masses and high mass-loss rates, and 
is supported by the evaporation time scales of the order of 
10'' yr computed by Henney & O'DeU (1999). One would 
expect that larger disks also lead to the formation of more 
extended proplyds. Larger disks could, e.g., be the result 
of an initially higher angular momentum in the collapsing 
molecular cloud core. 

With the aim to simulate the physical conditions close 
to those observed in NGC 3603 for Proplyd 1, the EUV 
photon rate was set to 10^^ photons s^^, and the distance 
to the ionizing source to 4.01 x 10^^ m (1.3 pc). The effec- 
tive temperature of the central source was set to 38,500K 
and the bolometric luminosity to 2.02xlO^Lo. 

We carried out two sets of simulations for different FUV 
photon rates. For the first set the FUV photon rate was de- 
termined from a blackbody spectrum. As indicated in Ta- 
ble ^, this leads to an underestimate of the true FUV pho- 
ton rate and results in lower mass-loss rates and a smaller 
size of the proplyd. For the second set of simulations the 
FUV photon rate was set to 1.2x10^^ photons s~^. In the 
following we will discuss only the results from the second 
set of simulations. The initial conditions for this set of 
simulations are summarized in Table 0. 

5.2. Mass-loss rates and life expectancy 

Figure |^ shows the evolution of the proplyd with time 
from the instant the Lyman continuum flux was turned on. 
An ionization front engulfing the star-disk-envelope system 
develops almost instantaneously. Initial mass-loss rates of 
the order of 10~^ Mq and evaporation flow velocities of the 
order of 20kms~^ lead to a rapid depletion of the central 
mass reservoir. After 50,000 yr the disk-envelope system 

^Note that the disk radii r,^ given by Johnstone et al. (1998) arc 
in the column heads of Table 1 and 2 in their paper. 



has already lost 2/3 of its inital mass. After 100,000 yr 
only 0.1 M0 remains in the disk. After 140,000yr almost 
95% of the disk mass has been evaporated. 

The simulations confirm the finding by Johnstone et al. 
(1998) that the FUV photon rate drives the mass-loss by 
heating up the region between the embedded disk and the 
ionization front. The resulting neutral flow influences the 
size of the ionized envelope. The simulations agree well 
with time scale estimates for disk photoevaporation by 
HoUenbach, Yorke, & Johnstone (1999). 

5.3. Emission-line maps in Ha 

In an attempt to compare the model to our HST obser- 
vations, we computed the emission line maps in Ha. The 
map shown in Figure ^ is based on the numerical model at 
a time of 50,000 yr. A viewing angle of 60° was assumed. 
The map was convolved with a theoretical HST/WFPC2 
point spread function computed with TinyTim (Krist & 
Hook 1997) and then resampled to the resolution of the 
Planetary Camera. 

Overall, the simulated Ha map shows a good resem- 
blance to the observations, and the models provide a viable 
explanation for the proplyds in NGC 3603. 

Despite the high FUV photon rate, the resulting pro- 
plyd is still a factor of two smaller than the proplyd-like 
structures observed in NGC 3603. This would indicate 
that either the central disks are larger, or that the pro- 
plyds in NGC 3603 have only very recently (within the 
last lO'^yr) been exposed to the UV radiation field. The 
latter would cause timescale problems, as it would require 
a very rapidly receeding HI/HII front. 

6. SUMMARY AND OUTLOOK 

While the present analysis cannot give a definite answer 
about the nature of the proplyd-like nebulae in NGC 3603, 
the WFPC2 data already provide a number of important 
clues: 

i) The tadpole shape suggests that the dynamics of the 
proplyds are dominated by the central cluster of NGC 
3603. 

ii) The variation of the Ha brightness of the proplyds 
with distance to the ionizing cluster further indicates that 

actually in units of 10"'^'' cm, not in units of 10^^ cm as erroneously quoted 
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Fig. 4. — Comparison of an Ha emission line map of the 50,000 yr old model and observed proplyds. The model has been 
convolved with a TinyTim Point Spread Function (Krist & Hook 1997) and resampled to the pixelscale of the Planetary 
Camera. 



the proplyds are ionization bounded. 

iii) The HST/WFPC2 and VLT/ISAAC images nicely 
resolve the nebulae into a region of neutral material 
surrounded by an outer ionization front. As discussed 
by O'Dell (1998), the fact that proplyds are ionization 
bounded implies that no EUV photons, but only FUV 
photons penetrate and heat up the central region of the 
proplyds. The neutral interior provides a steady supply of 
material for the ionization front. 

iv) The key to a better understanding of the internal 
density structure of the nebulae is the direct detection of 
central point sources, such as the source possibly associ- 
ated with Proplyd 1. The near- infrared photometry alone, 
however, does not really allow for a good estimate of the 
overall luminosity of the embedded source and the amount 
of internal extintion. Similarly, the limited spectral reso- 
lution of our ground-based spectra does not allow us to 
determine mass-loss rates for the proplyds. 

v) Faint nebular arcs detected between the proplyds and 
the cluster can be explained by bow shocks resulting from 
the interaction of the evaporation flow with the winds from 
the high-mass stars in the cluster. 

The simulations indicate that high mass-loss rates of the 
order of 10~^ Mq yr~^ are necessary in order to explain 
the physical size of the proplyds. Similar to the Trapez- 
ium system, the proplyd-like structures in NGC 3603 ap- 
pear to be short-lived phenomena. The fact that we de- 
tect three proplyds in NGC 3603 then suggests that pro- 



plyds/circumstellar disks are common not only in Orion, 
but throughout the Milky Way. 

Another consequence of the short life expectancy of pro- 
plyds in giant HH region is that they are not likely to 
survive long enough to form planetary systems, which are 
expected to form on timescales of the order of 10^ yr (Lis- 
sauer 1987). It might, however, also be possible to form 
giant planets on a considerably shorter timescale (Boss 
1998). 

While we do not have definite proof that the proplyd-like 
features in NGC 3603 are of the same nature as the pro- 
plyds in Orion, the similarities in morphology and phys- 
ical characteristics strongly suggest that they are related 
phenomena. The numerical modelling shows that the pro- 
plyds can indeed be scaled to the size and the physical 
conditions in NGC 3603. 

Further observations are needed to elucidate the nature 
of the proplyd-like nebulae in NGC 3603. High-spatial res- 
olution imaging and spectroscopy with HST in the optical 
and near-infrared will provide additional insights. Similar 
to Orion, it might be possible to detect the central disks 
glowing in the optical [Oi] emission line, and in the near- 
infrared molecular hydrogen lines. The physical extent, 
location within the proplyd, and inclination of the disk 
would provide important constraints for the simulations. 

While the radio continuum survey by De Pree et al. 
(1999) at 3.4 cm did not detect the proplyds, in the future 
it might be possible to detect and resolve the disks with 
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the Atacama Large Millimeter Array (ALMA), which is 
currently under study for the Chilean Atacama desert. 

Deep ground-based thermal infrared imaging with the 
new generation of 6 to 8 m class telescopes in the south- 
ern hemisphere could lead to the detection of more heavily 
embedded protostars in the center of the proplyds. The 
luminosity of each protostar would provide important con- 
straints on its mass and evolutionary status. 
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Table 1 

Location and size of the three proplyd-like structures 

Name RA(2000) DEC(2000) distance from cluster 



size 



Proplyd 1 


Ilhl5ml3.13s 


-61° 15' 50.0" 


43.6" (1.3 pc) 


1.8"" X 3.2" 


Proplyd 2 


llh 15m 16.59s 


-61° 16' 06.2" 


72.5" (2.2 pc) 


1.4" X 3.5" 


Proplyd 3 


llh 15m 07.73s 


-61° 15' 16.8" 


68.0" (2.0 pc) 


1.2" X 2.5" 



a:0.9", if we assume that Proplyd 1 is a superposition of two or three individual proplyds 
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Table 2 

Surface brightness of proplyd heads in Ha (F656N), [Nil] (F658N), V (F547M), R (F675W), and I (F814W) measured for 
a 0.5" aperture radius in the HST VEGAMAG system. 

Name Ihc [cm-^s"^]" Ha [mag/"^] [Nil] [mag/"^] F547M [mag/"^] F675W [mag/"^] F814W [mag/"^] 



Proplyd 1 
Proplyd 2 
Proplyd 3 



0.56 
0.20 



13.60 
14.73 



15.90 
16.71 



24.84 



21.58 



20.91 



a: Flux corrected for an assumed foreground extinction in Ha of AHa=4mag 



Table 3 

Near-infrared surface brightness of proplyd heads and photometry of nearby point sources "*" 
Name J. H K, 



(VLT/ISAAC). 



Proplyd 1 
Proplyd 2 
Proplyd 3 



Proplyd 1* 
Proplyd 3* 



14.8 mag/"^ 
16.2 mag/"2 
14.0 mag/"^ 



13.7 mag/"^ 

15.7 mag/"2 

12.8 mag/"^ 



15.4 mag 
12.9 mag 



14.1 mag 

13.2 mag 



12.9 mag/"^ 
14.5 mag/"2 
12.2 mag/"^ 



13.4 mag 
12.6 mag 



Ratio of FUV to EUV photon rates 
Toff [K] blackbody stellar photosphere" 
Z = O.2Z0 Z = Zq 



30,000 


4.9 


45.5 


130 


35,000 


3.1 


6.4 


6.9 


40,000 


2.2 


2.3 


2.4 


45,000 


1.6 


1.7 


1.4 



Table 4 



a: 



Derived from stellar photosphere models by A. Pauldrach et al. (1998). 



Table 5 



Model Parameters 

distance to ionizing source 

Tcff 

Luminosity 

EUV flux (hi/ > 13.6 eV) 

FUV flux (6eV < hi^ < 13.6 eV) 



4.01 X 10^6 m (1.3pc) 
38,500K 

2.02 X 10^ Lq 
10^1 s-i 

1.2 X 10^2 g-i 



